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ABSTRACT 

We study the formation histories and present-day structure of satellite galaxies formed 
in a high resolution hydrodynamic simulation of a Milky Way-like galaxy. The simu- 
lated satellites span nearly 4 orders of magnitude in luminosity but have a very similar 
mass within their inner 600 pc, ~3x 10 7 Mq, with very little scatter. This result is in 
agreement with the recent measurements for dwarf spheroidal galaxies (dSphs) in the 
Milky Way by Strigari et al. In our simulations a preferred mass scale arises naturally 
from the effects of the early reionisation of gas. These impose a sharp threshold of 
~ 12 km s _1 on the circular velocity of haloes which can cool gas and make stars. 
At the present day, subhaloes that host satellites as luminous as the classical Milky 
Way dwarfs (Ly ^ 2.6 x 10 5 Lq), have typically grown to have circular velocities 
of > 20 km s _1 . There are, however, subhaloes with similar circular velocities today 
which were, nevertheless, below threshold at reionisation and thus remain dark. Star 
formation in above-threshold haloes is truncated when the halo is accreted into the 
main galaxy progenitor. Thus, many properties of today's dwarf satellites such as their 
luminosity and star formation rate are set by their accretion time. 

Key words: methods: numerical - galaxies: evolution - galaxies: formation - cos- 
mology: theory - galaxies: dwarf. 



1 INTRODUCTION 

The mass density of the Universe is dominated by exotic 
dark matter whose identity remains unknown. Important 
clues could be hidden in the structure of the dwarf satellite 
galaxies of the Milky Way which are the most dark matter 
dominated systems known in the Universe. It has recently 
been found that these galaxies all seem to have approxi- 
mately the same central densit y (with i n 300 or 600 pc), ir- 
respective of their luminosity (|Mateol [1991 iGilmore et al.l 
120071 ; IStrigari et alj|2007l , |2008| ). There thus appears to be a 
special scale or threshold involved in the formation of dwarf 
satellites. A threshold central density could arise naturally 
if there is a cutoff in the primordial power spectrum of den- 
sity fluctuations produced by free-st reaming of warm dar k 
matter particles in the early universe IjBond fc Sza lav 1983). 
However, if the dark matter is cold, then there are no pri- 
mordial processes that operate on the relevant scales. In this 
case the origin of the threshold must be a late-time astro- 
physical process that inhibits star formation in small dark 
matt er haloes, such as the reionisation of gas at early t imes 
(e.g. lEfstathiou|[l992l : lOkamoto. Gao. fc Theunsll2008al ) . 
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By suppressing star formation in small cold dark mat- 
ter haloes, early re ionisation would al s o solve the 'missin g 
satellite' problem jKlvpin et alj Il999l : iMoore et al.lll999T l. 
as demo nstrated using semi-analytic models of galaxy for- 
mation dBullock et al lhood ; [Benson et al.ll2002l:lSomervillel 
120021: iKravtsov et alj|2004l ). Recently. iMaccio et alj (|2009n 
and I Li et al. I (|2009l) have confirmed, also with semi-analytic 
techniques, that this process sets the common mass scale 
seen in the inner 300 (or 600) pc of dwarf sat ellite 
haloes. One problem with the work of iLi et all (|2009l ) and 
M accio et alll|2009l) . however, is that they modelled reion- 
isation using the filtering mass approach of iGnedinl (2000) 
which has been shown to overestimate the impor tance of 
this effect (|Hoeft et al.ll2006l ; lOkamoto et al.ll2008al ). 

Since the reionisation of gas at early times and its ef- 
fects on star formation are complex processes, full hydrody- 
namic simulations are required to determine how they influ- 
ence the structure and evolution of dwarf galaxies. In this 
Letter, we analyse a high resolution hydrodynamic simula- 
tion that follows the formation of a Milky Way-type galaxy 
and is well suited to a detailed study of the formation of 
dwarf satellite galaxie s. This simulation has been shown by 
Oka moto et al. I (|2009l . Fig. 10) to give an excellent match 
to the observed number of faint satellites in the Milky Way 
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down to the resolution limit of the simulation (V-band mag- 
nitude Mv = —7), although it predicts somewhat too few 
very bright {Mv < —18) satellites. Here, we focus on the 
physical processes that determine whether or not a subhalo 
is able to form a galaxy and on the inner structure of the 
resulting object. 



2 THE SIMULATIONS 

We study the satellite population that formed in a high- 
resolution resimulation, using full baryonic physics, of 
'Aquarius-D', one of the 6 cold dark matter haloes simulate d 
as part of the 'Aquarius' project bv ISpringel et all (2008). 
Our simulation includes all the gas physics thought to be 
relevant for galaxy formation: a time -evolving photoionising 
background ( Haardt fc M adau 200 ]]), metallicity-depe ndent 



gas coolin g and photoheating JWiersma et al.l l2009h . SN 
feedback (Oka moto et al.l l2008bl ). and chemical evolution 



due to SNe II, SNe la, and AGB stars ilOkamoto et al.ll2005h. 
Detail s of this simulation may be found in Ok amoto et al.l 
(|2009l . their halo 'Aq-D-HR'). 

For our analysis two processes are particularly impor- 
tant: reionisatio n and galactic winds. T he former is treated 
as described in lOkamoto et aL ifeOOSa l '). Fo r the latter, we 
use the model introduced bv lOkamoto et al.l (|2009h . In brief, 
when a gas particle receives an amount of energy AE dur- 
ing a time-step At, this particle is added to a wind with 
probability, p w = A_Z5/(im g asUw), where u w is the initial 
wind speed. This is given as v w = 5<r, where a is the 
one-dimensional velocity dispersion of dark matter parti- 
cles around the gas particle. This expression implies that 
the wind mass generated by an SN is proportional to a~ 2 , 
i.e. per unit of star formation, less massive galaxies blow 
more (but slower) winds. This model reproduces both the 
observed luminosity function of the Milky Way satellites and 
the luminosity-metallicity relation. 

The mass of dark matter and gas particles in the sim- 
ulation are ttidm — 1-5 x 10 5 H~ 1 Mq and m gas ~ 3.5 x 
10 4 ft _1 M0, respectively. The gravitational softening length 
for the dark matter, gas, and star particles is kept fixed in 
comoving coordinates for z > 3; thereafter it is frozen in 
physical units, e = 0.175 h~ 1 kpc. To test for numerical con- 
vergence, we also ran a low resolution version of the simula- 
tion with m D M ^ 1.9xl0 6 /i _1 M , m gas ~ 4.1xl0 5 h~ 1 M e , 
and e = 0.425 /i _1 kpc. 

In order to identify subhaloes and satell ites, we use 
the SUBFIND algorithm (|Springel et alj|200ll ). We define 
subhaloes as systems that consist of at least 32 particles, 
and satellites as subhaloes that contain at least 10 star 
particles, unless otherwise stated. Satellite luminosities are 
computed using the population synthesis code PEGASE2 
jFioc fc Rocca-Volmerangd[l997j '). 



3 RESULTS 

The main result of our study is presented in Fig. [T]where, in 
the upper panel, we plot the mass within the inner 600 pc, 
Mo6, of the 38 satellites resolved in the simulation, as a func- 
tion of the IZ-band luminosity. The mass varies by less than 
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Figure 1. Upper panel: mass within the inner 600 pc of satellite 
galaxies as a function of V-band luminosity. The stars and squares 
show the mass obtained directly from the high- and low-resolution 
simulations, respectively. The points with error b ars show mea- 
surem ents for the classical dSphs of the Milky Way l lStrieari et al.l 
2007). The values of Mga derived from the empirical relation be- 
tween Mo6 and v Cj o are indicated by the plus signs. Lower panel: 
as the upper panel, but for the maximum of the circular velocity 
profile of satellites, « c ,o- 



a factor of 5 over 4 orders of magnitude in lum inosity, repro- 
ducing very well the overall behaviour of the IStrigari et al.l 
|2007ft data for classical dSphs, shown by the black dots. 
The lower panel shows the maximum of the circular veloc- 
ity curve, v C: o = u c (r- m a X ), which also agrees well with the 
data. The latter agreement is perhaps not too surprising 
given that in order to deriv e value s of -u c ,o for the Milky 
Way satellites, IStrigari et all (|2007h adopted a relation be- 
tween v c ,q and r max de rived from the ACDM simulation of 
iDiemand et~aD (|2007bh . 

Before examining the physical processes responsible for 
the result of Fig.[T] we discuss the reliability of our numerical 
data. Both panels of this figure include results from our low 
resolution simulation. For u Cj o these agree very well with 
those from the high resolution simulation, indicating that 
this quantity has converged. For Moe on the other hand, the 
results from the low resolution simulation are about a fac- 
tor of 2 lower than for the high resolution simulation. The 
explanation of this apparent paradox is simply that, v c p 
is numerically more robust than Moe (ISpringel et alj|2008l ). 
To test whether or not Moe has converged in our high res- 
olution simulation, we make use of the empirical relation 
between Moe and v c ,o derived from the highest resolution 
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Figure 2. The cumulative distribution of maximum circular ve- 
locity for all subhaloes and visible satellites. The solid line gives 
the number of dark matter subhaloes in the halo and the solid cir- 
cles show the number of simulated satellites that are brighter than 
Draco. The distribution for the full satellite sample (iV s tar 10) 
is shown with open circles. The black arrow indicates the value of 
the circular velocity (23.5 km s —1 ) below which failed subhaloes, 
i.e. those that contain no star particles in the simulation, exist. 

Aquarius simulatiorQ (which has about 100 times better res- 
olution than o u r high resolution hydrodynamic simulation; 
ISpringel et ail (|2008T )). As shown in Fig. [TJ the values of 
Mo6 derived from this relation agree well with the directly 
measured values, showing that the high resolution simula- 
tion gives converged and reliable values of Mo6 for all the 
satellites resolved in the simulation. 

The lower panel of Fig. fTJ provides an important clue 
to the reasons behind the preferred mass scale picked out 
in the top panel. Over the range spanned by the observed 
dSphs (Lv < 2 x lO 7 Z/0), the velocity u Ci o is essentially 
constant, at ~ 20 km s~ , and then it rises slowly with 
luminosity at brighter magnitudes. Thus, there appears to 
be a threshold circular velocity required for a dwarf satellite 
to form. This threshold is evident in Fig. [2] which compares 
the circular velocity distribution of satellites to that of the 
entire subhalo population. All subhaloes with v c .o > 25 km 
s _1 succeed in making a satellite, but below this value, the 
fraction of successful subhaloes drops rapidly. This inability 
of small subhaloes to form a visible galaxy is the reason 
why there is no 'satellite problem' in the CDM model: only 
a small fraction of the large subhalo population are lit up 
by galaxy formation. 

Note, however, that there is no hard cutoff in v c .o in 
Fig. [2] Below f c ,o = 23.5 km s~ , there are subhaloes that 
succeed in making satellites as luminous as the classical 
dSphs (Lv ^ 2.6 x 10 5 Lq), but also failed subhaloes that 
do not contain any stars. The difference between haloes of 

1 logKo/(km s- 1 )] = 0.581og[Moe/(/i -1 M©)] - 3.092 
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Figure 3. The maximum circular velocity, n c ,g, at the epoch of 
reionisation, z rc = 9, of the progenitors of present-day subhaloes, 
plotted against their maximum circular velocity today, v Cj o- Satel- 
lite galaxies containing at least ten star particles at z = are 
represented by the solid circles while dark satellites that do not 
contain any star particles at z = are shown by squares. We also 
show subhaloes containing at least two star particles at z = 
by open circles. There is a threshold value of v c ^q that separates 
successful from failed subhaloes. For reference, we represent the 
threshold by the dotted line at v C: g = 12 km s . 

similar, relatively large circular velocities which do and do 
not form a visible galaxies is apparent in Fig. [3] Here we 
compare the present day maximum circular velocity of a sub- 
halo, v c ,o, with the maximum circular velocity of its main 
progenitor, u c ,9, at the epoch of reionisation, which in our 
simulations happens at z la — 9. This figure shows a sharp 
transition in « c ,g between subhaloes whose descendants go 
on to make visible galaxies and those whose descendants do 
not. Surprisingly, this threshold applies not only to lumi- 
nous satellites, but also to subhaloes that contain only two 
star particles. The transition occurs at v c ,9 — 12 km s _1 and 
is due to the effects of reionisation: the increase in thermal 
pressure which follows re ionisation evaporates the gas from 
subthreshold haloes (e.g. lReeslll986l ; lOkamoto et alj|2008ah . 

Haloes with maximum circular velocity above the tran- 
sition critical velocity, v Ct g > 12 km s _1 , are observed as 
dwarf satellites today. Their circular velocities today peak at 
around « Cj o ~ 20 km s _1 , as may be seen in Fig. [3] This is the 
origin of the common velocity scale seen in the lower panel 
of Fig. [1] This characteristic velocity translates directly into 
the characteristic inner halo mass mea s ured by observers in 
today 's satellites ()Gilmore et al ] |2007l : IStrigari et al 

boost) . 

In our simulations, reionisation is the main process that 
separates successful from unsuccessful haloes. However, the 
properties of a satellite are determined by subsequent events. 
Particularly important is the accretion of the (sub)halo into 
the main progenitor of the host galaxy. The epoch of ac- 
cretion determines many of the key properties of the final 
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(a) dwarf satellites 




(b) dark satellites 
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Figure 4. Formation histories of dwarf satellites and dark sub- 
haloes, (a) total mass (solid), gas mass (dotted) and stellar mass 
(dashed) in the main progenitors of a selection of simulated haloes 
that end up as dwarf satellites today, plotted as a function of 1 + 2. 
Reionisation occurs at 1 + z Ie = 10. Arrows indicate the redshift 
when the progenitor is accreted into a larger halo (i.e. when it 
becomes a subhalo). Lower panels: the maximum of the circu- 
lar velocity of the haloes. The horizontal dotted line marks v c> g 
while the horizontal dashed line indicates the critical velocity at 
the epoch of reionisation (v c ,g = 12 km s ). (b) as (a), but for 
failed subhaloes that end up as dark satellites today. 



satellite, such as its luminosity and star formation rate. The 
key events in the life of a subhalo are illustrated in Fig. [4] 

Consider first the evolution of haloes that go on to 
make visible satellites by the present day, illustrated in pan- 
els (a). As may be seen in the lower subpanels, at z rc — 9 all 
these haloes have above-threshold maximum circular veloc- 
ity (v c ,9 > 12 km s _1 ). The gas in these objects, shown by 
the red line in the upper subpanels, is unaffected by reion- 
isation and these haloes are able to form stars (blue lines). 
Eventually, however, the halo is accreted into the progenitor 
halo of the main galaxy, becoming a subhalo. At this time, 
the remaining gas is stripped completely and star formation 
ceases. This explains why most of today's dwarf satellites 
are old. Although we do not find the correlation between 
accretion redshift and halo mass at t he epoch of accretio n 
reported in the semi-analytic models of lMaccid et al.ll|2009h . 
we do find a correlation between accretion redshift and final 
satellite luminosity. The satellite illustrated in the 3rd col- 
umn has a slightly unusual behaviour. This object loses its 
gas well before becoming a subhalo. Since its maximum cir- 
cular velocity at z = 9 is just above the critical value, reion- 



isation and photoheating do heat up the gas which becomes 
loosely bound and is subsequently blown out by supernova 
explosions and increased background radiation by z ~ 6. 
The photoionising background prevents gas f rom reaccreting 
into a halo with such small circular velocity (|Okamoto et al.l 
l2008al ). 

Panel (b) follows the fate of typical failed satellites, i.e. 
of haloes that are below the circular velocity threshold at 
reionisation and are destined to remain as dark subhaloes 
today. Most of them lose all of their gas by evaporation at 
z Ic . Although they can subsequently gain small amounts of 
gas as they grow in mass (and circular velocity), the gas 
density in their shallow potential wells is much too low for 
stars to form. The halo illustrated in the rightmost panel has 
value of Uc,9 which is only just above threshold. As a result, 
it does not lose its gas immediately after reionisation but 
holds on to it for some time before finally losing it at z ~ 2. 
This object never manages to make stars because the large 
thermal pressure prevents gas from attaining a sufficiently 
large density and so the (sub)halo remains forever dark. 



4 SUMMARY AND DISCUSSION 

In this Letter we have used a high resolution N-body hy- 
drodynamic simulation to investigate the formation histo- 
ries of the satellites that form in a Milky- Way type galaxies 
in a ACDM universe. These simulations have been shown 
to reproduce the main optical and chemical properties of 
the Local Group satellites , including their luminosity func- 
tion |Okamoto et al. 2009). In this study we have focused on 
the properties that distinguish the small subset of subhaloes 
that succeed in making a galaxy, from the vast majority that 
remain dark. 

Our simulations naturally reproduce the observation 
that the Milky Way satellites have very similar central den- 
sities irrespective of their luminosity. They also show that 
there is a threshold circular velocity today, v Cl o — 23 km 
s~\ above which all subhaloes present today host satellites 
with the luminosities typical of the classical satellites of the 
Milky Way (L v > 2.6 x 10 5 L ). The number of satellites 
drops rapidly below this critical circular velocity but around 
this value there are haloes which succeed in making a galaxy 
and haloes which do not. The distinction between them can 
be traced back to the epoch of reionisation when there is a 
hard threshold in maximum circular velocity, v c ,9 — 12 km 
s _1 , that divides haloes where gas can cool and make stars 
from haloes in which it cannot. The scatter between v c ,o and 
v c ,9 that develops over time explains why seemingly similar 
subhaloes today can have such different mass-to-light ratios. 
This key feature cannot be explained by the effect of galac- 
tic winds alone. Galactic winds in which the wind velocity 
scales in proportion to velocity dispersion were shown by 
jOkamoto et al.ll2009T l to be required to explain the satellite 
luminosity function. Here we have shown that reionisation 
is required to explain the circular velocity function. 

Haloes with sub-threshold u Ci g lose their gas by evap- 
oration during reionisation while the haloes with above- 
threshold n c , 9 retain their gas which eventually overcomes 
the increased thermal pressure produced by reionisation and 
proceed to make stars. This explains why lOkamoto et al.l 



The origin of dark satellite galaxies 5 



(2009) could not find any residual signs of reionisation in 
the star formation histories of simulated satellites. 

If the CDM model is right, a large number of small, 
completely optically dark subhaloes should be orbiting in 
the Milky Way halo. They might one day be detected in- 
directly either through annihilation radiation if the dark 
matt e r consists of supersymmetric parti c les (e.g. | Baltz et all 



200d; ICalcaneo-Roldan fc Moord 120001; lEvans et al.1 1200I 



Diemand et al.ll2007al ; ISpringel eta l. 2008) or through grav - 



itational lensing effects (e.g |Kochanek|[l99ll ; IXu et all 2009) 
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Broadly speaking, our results agree with those obtained 
using semi-analytic models applied to high-resolution N- 
body simulations. These studies also find that a preferred 
mass scale arises naturally in dwarf galaxies from a com- 
bination of astrophysical processes that inhibit star forma- 
tion in small haloes and the fact that the mass within 600 
pc (or 300 pc) is not a very sensitive function o f the total 
subhalo mass (|Li et alj 12009; iMaccio et alJ l2009L It is nev- 
ertheless surprising that lLi et aTT feoog) managed to obtain 
a correlation between Moe and satellite luminosity similar 
to ours (Fig. [TJ and to the real data because they assumed 
an extreme reionisation model w hich is incompatible w ith 
the results of recent simulations l|Okamoto et al.l|2008al ). It 
is possible that the minimum mass scale in their model was 
set by their assumption that gas does not cool in haloes 
with virial temperature below 10 4 K rather than by the ef- 
fects of photoionising background. Note that the threshold 
circular velocity imposed by reionisation in our simulations, 
v c ,9 — 12 km s~ , corresponds to a virial temperature of 
~ 5 x 10 3 K. 

Using a simple analyti cal model applied to an N-body 
simulation. iKravtsovl (|2009l ) also finds that a preferred satel- 
lite mass scale (as well as other properties of the satellite 
system, such as the luminosity function and spatial distri- 
bution) can be explained if the star formation efficiency 
decreases with decreasing halo mass. In our simulations a 
dependence of star formation efficiency on halo mass of this 
kind is a natural outcome of the behaviour of galactic winds. 
Kravtsov's model does not assume a sharp threshold in the 
star formation efficiency. This, however, arises naturally in 
our simulations as a result of reionisation and seems to be 
required in order to explain why some dark matter halos 
always remain dark. 

Our simulations ignore the possibilities that gas may 
be self-shielded against ionising radiation, or that star for- 
mation may result from molecular gas cooling before reion- 
isation. Both of these proc esses could be important in the 
formation of small objects (ISusa fe Umemu ra 2004a, bj) and 
may allow small amounts of star formation in sub-threshold 
haloes. It will be particularly interesting to explore how im- 
portant t hese effects are for th e newly-discovered ultrafaint 
satellites l|Koposov et al.ll2008l ) . We leave these issues for fu- 
ture work. Our present simulations, however, demonstrate 
that a preferred mass scale in satellites does not require 
abandoning the standard model of cosmology but is, in fact, 
a natural byproduct of galaxy formation in small cold dark 
matter haloes. 
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